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Abstract* 

The  problem  of  selecting  a  subset  containing  all  populations 
better  than  a  control  under  an  ordering  prior  is  considered.  Three 
new  selection  procedures  which  satisfy  a  desirable  basic  requirement 
on  the  probability  of  a  correct  selection  are  proposed  and  studied. 
Two  of  the  three  procedures  use  the  isotonic  regression  over  the 
sample  means  of  the  k-treatments  with  respect  to  the  given  ordering 
prior.  Tables  of  constants  which  are  necessary  to  carry  out  the 
selection  procedures  with  isotonic  approach  for  the  selection  of 
unknown  means  of  normal  populations  are  given.  The  results  including 
Monte  Carlo  studies  indicate  that,  in  general,  the  stepwise  procedure 
based  on  isotonic  estimators  is  the  best. 
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1 .  Introduction 

In  this  paper,  three  new  selection  procedures  are  given  for  the 
problem  of  selecting  a  subset  which  contains  all  populations  better 
than  a  standard  or  control  under  simple  or  partial  ordering  prior. 

Here  by  simple  or  partial  ordering  prior  we  mean  that  there  exist 
known  simple  or  partial  order  relationships  (defined  more  specifically 
later  in  Section  2)  among  unknown  parameters.  The  procedures  described 
do  meet  the  usual  requirement  that  the  probability  of  a  correct 
selection  is  greater  than  or  equal  to  a  predetermined  number  P*,  the 
so-called  P*-condi tion. 

Many  authors  have  considered  the  problem  of  comparing  populations 
with  .1  control  under  different  types  of  formulations  (see  Gupta  and 
Paf.chapakesan  (1979)).  Dunnett  (1955)  considered  the  problem  of  sepa¬ 
rating  those  treatments  which  are  better  than  the  control  from  those 
that  are  worse.  Gupta  and  Sobel  (1958),  Gupta  (1965),  Nalk  (1975), 

Brostrom  (1977)  studied  the  problem  of  selecting  a  subset  containing 
all  populations  better  than  the  control.  Lehmann  (1961)  discussed 
similar  problems  with  emphasis  on  the  derivation  of  a  restricted  mlnlmax 

procedure.  Gupta  and  Kim  (1980),  Gupta  and  Hsiao  (1980)  studied  the  problem  of 

*This  research  was  supported  by  the  Office  of  Naval  Research  contract 
N00014-75-C-0455  at  Purdue  University.  Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of  the  United  States  Government. 
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selecting  populations  close  to  a  control.  In  all  these  papers  It  Is 
assumed  that  all  populations  are  Independent  and  that  there  Is  no  In¬ 
formation  about  the  ordering  of  unknown  parameters.  However,  In  many  sit¬ 
uations,  we  may  know  something  about  the  unknown  parameters.  What  we 
know  Is  always  not  the  prior  distributions  but  some  partial  or  incom¬ 
plete  prior  Information,  such  as  the  simple  or  partial  order  relation¬ 
ship  among  the  unknown  parameters.  This  type  of  information  about  the 
ordering  prior  may  come  from  the  past  experiences;  or  It  may  arise  in 
the  experiments  where,  for  example,  higher  dose  level  of  a  drug 
always  has  larger  effect  on  the  patients. 

In  Section  2  definitions  and  notations  used  In  this  paper  are 
introduced.  In  Section  3  we  consider  the  problem  for  location  param¬ 
eters.  We  propose  three  types  of  selection  procedures  for  the  cases 
when  the  control  parameter  is  known  or  not  known  (the  scale  parameter 
may  or  may  not  be  assumed  known).  Some  equivalent  forms  of  the  pro¬ 
cedures  are  given,  and  their  properties  are  discussed.  In  Section  3 

simple  ordering  priors  are  assumed  and  some  theorems  in  the  theory 
of  random  walks  are  used.  A  selection  procedure  for  the  problem  of 
selecting  all  populations  better  than  the  control  under  partial  ordering 
prior  is  given  in  Section  4.  Section  5  deals  with  the  use  of  Monte  Carlo 
techniques  to  make  comparisons  among  the  selection  procedures  proposed 
in  Section  3  and  those  in  Section  4,  respectively. 
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2.  Notations  and  Definitions 

Suppose  we  have  k  +  1  populations  ttq,  The  population 

treatment  tt0  Is  called  the  control  or  standard  population.  Assume 
that  the  random  variable  X^is  associated  with  F( •  ;0^ )  and  X^,...,X^n  , 

i  =  l,...,k,  are  independent  samples  from  it^  Assume  that  we  have  an 

ordering  prior  of  @^,...,0^.  First  we  assume  that  the  ordering  prior 
is  the  simple  order,  so  that  without  loss  of  generality,  we  may  assume 
that,  0-j  .  £0^*  In  Section  4  we  will  consider  the  partial  order¬ 

ing  prior  case.  Note  that  the  values  of  e^'s  are  unknown. 

Suppose  our  goal  is  to  select  a  non-trivial  (small  )  subset  which  con¬ 
tains  all  populations  with  parameters  larger  (smaller)  than  the  control 
9q  (known  or  unknown)  with  probability  not  less  than  a  given  value  P*. 

The  action  space  G  is  the  class  of  all  subsets  of  the  set  {l,2,...,k}. 
An  action  A  is  the  selection  of  some  subset  of  the  k  populations,  i  6 A 
means  that  tt.  is  included  in  the  selected  subset. 

Let  £  =  (0q,  e^,...,e^).  Then  the  parameter  space  is  denoted  by  n, 

where  n  =  {£€F  j  e1  £  e2  £  . . .  <  e^;  -  ~  <  eQ  <  is  a  subset  of 

k+1 

k  +  1  dimensional  Euclidean  space  F 

The  sample  space  is  denoted  by  X  where 


X  =  (x  6  F 


nl+...+nk 


x  =  (x , X,  ,  . .  .  ,X.  ,  ,.  .  .X.  )}. 


(Here  0q  is  assumed  to  be  known). 

Definition  2.1.  A  (non-randomized)  selection  procedure  (rule)  6(x) 


is  a  mapping  from  X  to  G. 


A  population  ^  (1  =  1,...,k)  Is  called  a  good  population  If 
8^  >_  eQ.  A  correct  selection  (CS)  Is  the  selection  of  a  subset  which 
contains  all  good  populations.  A  selection  procedure  6  satisfies  the 
P*-condition  if 


Inf  Pfi(CS|6)  >  P*.  (2.1) 

c<£n  - 

Let  S  =  {6 1  Inf  P_  ( CS  I  <S )  >  P*}  be  the  collection  of  all  selection 

een  ! 

procedures  satisfying  the  P*-condit1on. 

In  the  sequel  we  will  use  the  isotonic  estimators  (see  Barlow, 
Bartholomew,  Bremner  and  Brunk  (1972)).  Hence  we  give  the  following  def¬ 
initions  and  theorems. 

Definition  2.2.  Let  the  set  7  be  a  finite  set.  A  binary  relation 
"<"  on  J  is  called  a  simple  order  if  it  Is 

(1)  reflexive:  x  <  x  for  x  6 J 

(2)  transitive:  x,  y,  z€J  and  x  <  y,  y  <  z  imply  x  <  z 

(3)  antisymmetric:  x,  y  and  x  <  y,  y  <  x  imply  x  *  y 

(4)  every  two  elements  are  comparable:  x,  y  6J"  Imply  either 

x  <  y  or  y  <  x. 

Definition  2.3.  A  partial  order  on  J  is  a  binary  relation  "<"  on  j,  such  that 
it  is  (1)  reflexive,  (2)  transitive,  and  (3)  antisymmetric.  Thus  every 
simple  order  is  a  partial  order.  We  use  poset  (Jt<)  to  denote  the  set 
J  that  has  a  partial  order  binary  relation  "<"  on  it. 
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Definition  2.4.  A  real-valued  function  f  is  called  isotonic  on  poset 
(J,<)  If  and  only  if  (1)  f  is  defined  on  j%  (2)  if  x,  y  €7,  x  <  y  imply 
f(x)  <  f (y) . 

Definition  2.5.  Let  g  be  a  real -valued  function  on  J  and  let  W  be  a 
given  positive  function  on  J.  A  function  g*  on  J  is  called  an  isotonic 
regression  of  g  with  weights  W  if  and  only  if: 

(1)  g*  is  an  Isotonic  function  on  poset  (J><) 

(2)  l  [g(x)  -  g*(x)]2W(x)  =  min  l  [g(x)  -  f(x)]2W(x), 

x€J  f€»  xfcJ 

where  3  is  the  class  of  all  isotonic  functions  on  poset  (J,<). 

From  Barlow,  et.  al .  (1972),  (see  their  Theorems  1.3,  1.6  and  the 
corollary  there),  we  have  the  following  theorems. 

Theorem  2.1.  There  exists  one  and  only  one  isotonic  regression  g* 
of  g  with  weight  W  on  poset  (J,<). 

There  are  some  known  algorithms,  such  as  the  "pool-adjacent-violators" 
algorithm  (see  page  13  of  Barlow,  et.  al.  (1972))  or  Aver,  Brunk,  Ewing, 
Reid  and  Silverman  (1955)  or  the  "up-and-dcwn  blocks"  algorithm,  Kruskal 
(1964),  which  show  how  to  calculate  the  isotonic  regression  under  simple 
order. 

The  following  max-min  formulas  were  given  by  Ayer  et.  al.  (1955). 

Theorem  2.2.  (max-min  formulas) 

Assume  that  we  have  poset  (J,<)  where  J  =  {0-j ,. . .  ,0^1,  0-|  <  . .  .<  0^, 
and  that  function  g:  3  -*■  F,  then  the  isotonic  regression  g*  of  g  with 
weight  W  has  the  following  formulas: 
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g*(oJ  *  max  min  Av(s,t) 
s<i  t>1 

*  max  min  Av(s,t) 
s<1  t>s 

=  min  max  Av(s.t) 
t>1  s<1 

=  min  max  Av(s,t) 
t>i  s<t 


where 


Av(s,t) 


l  g(e_)W(e  ) 
r=s  r  r 
T - 


I  H(er) 

r=s 


Corollary  2.1.  (g  +  c)*  =  g*  +  c,  (ag)*  =  ag*,  if  a  >  0,  c  €  IR  . 

Corollary  2.2.  [p(g*)g  +  ;p(g*)3*  =  p(9*)  +  ?(g*).  where  p  is  a 
nonnegative  function  and  cp  is  an  arbitrary  function. 


3 .  Proposed  Selection  Procedures  for  the  Normal  Means  Problem 

We  are  interested  in  the  (subset)  selection  problem  of  the  unknown 
means  of  k  normal  populations  in  comparison  with  a  standard  or  control 
normal  with  its  mean  known  or  unknown.  Thus  observations  are  taken  on 
which  are  independently  distributed  normal  random  variables 
N(n-,o  ),  j  =  l,...,n^;  i  =  l,...,k.  The  values  of  y i »y2  *  * ' *  are 

unknown,  but  their  ordering,  say,  '  -  yk  ^own.  Note  that 

in  this  case  we  replace  e  in  the  parameter  space  n  by  y,  all  other 
quantities  remaining  the  same. 


Let  us  define  the  subspace  a  (jj. €  o |  wk_.|  <  uQ  £  wk.1+1}  for 
1  ■  1,...,k-l,  the  subspace  nk  ■  {y  6  ft|  hq  <  H|)»  and  the  subspace 

k 

n0  “  *lL€nl  ^|<  <  *  ^en  we  ^ave  n  *  U  n^.  Note  that  the  control 

Uq  could  be  known  or  unknown.  If  Is  unknown,  we  assume  that  the 
distribution  of  population  ttq  Is  N(wq>  o^)  end  we  take  Independent 
observations  Xq-|  ....  .Xg^from  ttq  and  the  sample  space  X  becomes 

{X€lRn0+“‘+nk|  X  -  (XQ1 . X0no*  Xn . Xln  ,...,Xk1 . Xk  )}.  Usl 

the  partition  {fiQ,...,nkl  of  parameter  space  n,  we  have 

Inf  P  (CS|«)  =  inf  {inf  P  (CS|«)}, 
jj€n  ±  l<i<k  11 

for  any  selection  procedure  5£jB.  Hence  the  P*-cond1tion  Is  equivalent 


inf  P  (CS|6)  >  P*.  for  i  *  1 . k. 

u€iii  & 

Note  that  inf  P  (CS | « )  =  1  for  any  selection  procedure  5  since  there 

u€n0  — 

exists  no  good  population  in  this  case. 

Let  X.  =  xi  be  the  observed  sample  mean  from  population 
i  =  1 ,. . . ,k.  Let  J  denote  the  set  {y^ ,  u2>. . . ,yk>  where  w1  £  . . .  £  u^i 

and  let  W(y^)  =  n^o  =  w. ,  g(u^)  *  x^ ,  i  *  l,...,k.  Then  by  the  max- 
min  formulas,  the  isotonic  regression  of  g  is  g*.  where 

g*(y.)  =  max  min  “r - ,  1  =  l,...,k. 

l£s£i  s<t<k  p 

j  s  j 

The  isotonic  estimator  of  is  denoted  by  X^.k,  1  =  l,...,k  where 
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I 

f 

t 

i 


f 

i 


i 

I 

r 

f 


f 

t 

t 


i 

j 


r. 


i 


max  min 


t 

y  x.w. 
j»s  J  J 


lira*  III  I II  iT 

1  *K  l<s£i  s<t<k  |  w 
j-s 

»  max  {X...  } 


J 


1<s<1 


s:k 


(3.1) 


where 


*  ,  f„  xsws+xs+lwsfl  X£ws ■>•... +Xkwk 

xs:k  “  m1ntXs’  ws+ws+1  .  ws+...+wk  }>  (3,2) 

It  Is  known  that  the  Isotonic  estimators  fc.;k,  i  *  1,...,k  are  also  the 
maximum  likelihood  estimators  of  i  =  1,...,k. 

3.1.  Proposed  Selection  Procedure  6^ 

o 

Case  I.  u0  known,  common  variance  o  known,  and  common  sample  size  n. 
Definition  3.1.  We  define  the  procedure  6^  as  follows: 

Step  1.  Select  ir,,  1  *  l,...,k  and  stop,  if 

y  s  j(l)  ° 

Xl:k  -  u0  ’  dl:k 

otherwise  reject  and  go  to  Step  2. 

Step  2.  Select  ^ ,  i  *  2,...,k  and  stop,  if 

v  s  j(1  )  _o_ 

X2:k  i»0  '  d2:k  ^ 


otherwise  reject  itg  ancl  9°  t0  steP 


Step  k-1 . 


Select  ir.j ,  1  =  k-1,  k  and  stop,  if 


Xk-l:k  -  y0 


jd) 
Jk-1  :k 


ctep  K. 


otherwise  reject  rrk_^  and  go  to  Step  k. 
Select  *k  and  stop,  if 


otherwise  reject  irk> 

Here  dj]k's  are  the  smallest  values  such  that  6^  6  8,  that  is  6^  sat 
isfies  the  P*-condition. 


3.2.  On  the  Evaluation  of  Inf  P,  (C$|6,)  and  the  Values  of  the 
_ ktn,  11  1 

Constants  dj]),...,d£!/ 

For  any  u.€  1  £  i  <  k,  let  's  be  i . i .  d .  N ( 0 » 1 )  and  let  Zr.k 

m1n(Zr,  2  *»•••»  k-r+1  ' 

k"i+1  -  m  n 

■V  j=u, 

k-1+!  j  s  m  _ 

=  Pu(  U  U  (Xr.k  >M0  -  d Jv-SL}) 

^  j=l  r=l  r,k  0  Jn 


i  p  (  u  u  tzr;k  +  >  -  d \!J» 


j«l  r=l 


a/  /rf 


j:k 


which  is  increasing  in  y  ,  r  =  1 . k-1+1, 


Hence 


inf  Pj.tey,) 


On  the  other  hand, 


inf  P  (CS | 6, ) 
V  I 
u€n^  - 


k-1+1 


i  V  -  “0 


-  ■!'!  J" 


-0/7  >  _  ^  1 

PUk-i+l:k-  ak-i+l:k> 

whenever  y*  =  (yg,  -  -  ®,  Pg,.  . .  ,Pg)  €  n^. 


inf  P  (CS 1 61 )  =  P(Zk.1+1;k  1  -  dk!|+i:k)- 
y«^  - 


Thus,  we  have 
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Since  Z|<_^  +i  has  the  same  distributions  as  Z^. 
letting 


Vi 


(3.3) 


we  have 


It  is  clear  from  the  above  that  =  d|]|  for  all 

i  =  1,  2,...,k,  and  dj]j  is  increasing  in  i. 

p 

Theorem  3.1.  In  case  I,  (yQ  known,  common  known  a  and  common  sample 
size  n),  if  is  the  solution  of  equation 

P(Vi  >_  -  x)  =  P*  (3.5) 

where 

!  r 

V.  =  min  -  l  Z.  and  Z.  are  i.i.d.  N(0,1), 

1  l<r<i  r  j=l  J  1 

i  =  l,...,k, then  6^  satisfies  the  P*-condition. 

Proof.  For  any  i,  1  <_  i  <_  k, 

inf  P  (CS|5,)  =  P(V.  L  -  dPL.J  =  P*. 
y  1  1 

so  6-|  satisfies  the  P*-condition. 

Therefore,  the  problem  of  finding  the  djl^'s  reduces  to  finding  the 
distributions  of  V- |,...,V^.  This  is  achieved  by  using  some  results 
in  the  theory  of  random  walk. 
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3.3.  Some  Theorems  in  the  Theory  of  Random  Walk 

Suppose  Y.j ,  Y2>...  art  independent  random  variables  with  a  common 
distribution  H  not  concentrated  on  a  half-axis,  i.e.  0  <  P(Y-j  <  0), 
P(Y^  >  0)  <  1.  The  induced  random  walk  is  the  sequence  of  random 
variables 

S„=0,  S  =v,  *...+Yn>  ,=1,2,...  . 


(3.10) 


By  symmetry,  the  probabilities 

„n=P(S,  <0,....Sn<  0) 

have  the  generating  function  q  given  by 

log  q(s)  =  l  V  P^Sn  -  °^' 
n=l 

Note:  The  above  theorem  remains  valid  if  the  signs  >  and  <  are 
replaced  by  >  and  <,  respectively. 

Theorem  3.3.  The  generating  function  p(s)  of  P(Vj  >  x),  j  =  1,2,. 

is  1 

V  sj  P(V.  L  x)  =  exp  t  l  ~  sn  P(Sn>  0)1  (3.12) 

j=l  J  n=1 

where 

n 

Sn  =  £  U.  "  x),  n  =  1,  2 . 

n  1  =  1  1 

Proof.  Since  the  distribution  of  random  variable  Y.j  =  Z.  -  x  is  not 

r 

concentrated  on  a  half-axis,  and  Y^'s  are  i.i.d.  let  Sr  =  ^(Z^  - 

r  =  1 , . . . ,k.  Then 

{V.  >  x }  =  {  min  1  >  0}  =  {S,  >  0,...,S.  >01. 

0  "  r  r  1  J 

By  Teller's  Theorem  3.2  ,  we  complete  the  proof. 
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Now  let 


!.(*)  =  4,  -  P(S  »0)-  i  - 

J  J  J 


•<*>■  I  T  V 

n=l 


then  we  have 


p(s)  e  l  sj  P(V.  >  x)  -  exp  (a(s)). 
j=l  J 


_ ( n+1 ) f c )  =  ?  (*!)  p(j)(s)  a(n+1"j)(s),  for  all  n  >.  1 
Lemma  3.1.  P  V  H 

j=0 

Pro„f.  Since  p'(s)  -  P(.)  •  -suit  be  ^  ^  1n6,C 


tion  on  n. 


Theorem  3.4 


where 


Under  the  assumption  of  Theorem  3.3 


v  1  lim  _dn+  p(sj 
P(Vn+1  i  *>  *  Tn+TIT  sV  ds"+1 

/ 

•  i  I  P(VJ;>  x)Vj+t  "  =  0, 

n+  j=0  7 

p(Vn  >x)-=l,  for  a11  X- 


(3.13 


Proof.  By  Lemma  3.1,  we  have 


p(v  +1  i'x)  =  xhtttt  lim+ p(n  1)(s) 

s-*0 


*  lisWrrr^T! pU)(0)  Vl- 


1  ?  p!l1m  i  ^  . 

!i*tl  Tf 


■PTjlo  p(vi  -■  X)W 


14 


Let  6  (x)  -  P( V  >  x)  and  1-G  (x)  denote  the  limiting  distribution 
n  n  —  ® 

function  as  n  -*■  »  of  V  .  Suppose  the  distribution  of  random  variable 
-  x  is  not  concentrated  on  a  half  axis,  then  we  have  from 
Andersen-Fel ler  Theorem 

Gjx)  =  exp  {-  l  Jp(S  <0)1. 
r=l 

Now, let 

®J-  *  P*-  (3-'4) 


Now  we  can  use  the  recurrence  formula  of  Theorem  3.4  to 
solve  the  equations  P(V^  _>  -  djj  |+-|  *  p*,  1  = 


Remark  3.  1.  From  Section  3.2  we  know  that  dj^]+^.k  B  d]!|  0  *  l,...,k). 
The  values  of  dj|^,  for  k  =  1  (1)  6,  10,  «  and  P*  -  .99,  .975,  .95, 

.925,  .90  are  t’  ulated  in  Table  I. 

Definition  3.2.  We  define  a  selection  procedure  <$j  by  replacing 
the  inequality  in  the  ith  step  of  procedure  4.,  by  the  inequality 

V-k  -  "0  '  di:k  1  1,  ‘”k 

where  dj  •  •  >d£ ’k  are  sma^est  values  such  that  <sj  satisfies 
the  P*-condition. 


Then  it  can  easily  be  shown  that  the  selection  procedure  6^  and 
are  identical  and  dj]jj,  =  d!.^,  i  =  l,2,....k. 

3.4.  Some  Other  Proposed  Selection  Procedures  <52.  63,  <$4 
In  Case  I,  we  propose  some  other  selection  procedures: 


WJ  iiid  4HI  i  »-•  ,V  *-:  i  J 


hi  J  uU.'  l.i-1  »1 1  #1: 


Mi 
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Definition  3.3.  We  define  a  selection  procedure  fig  by 

fig:  Select  7^  if  and  only  if  X^  >_  mq  -  d  —  i  =  1 . .  ,k 

where  d  is  the  smallest  value  such  that  fig  satisfies  the  P*-condition. 

Note  that  under  assumptions  of  Case  I,  and  selection  procedure 
fig,  if  we  select  population  ,  then  we  will  select  populations 
for  all  j  ^  i,  since 


Evaluation  of  the  d-Values  of  fig 

For  any  i,  1  _<  i  <  k,  we  have  from  a  similar  argument  as  for  fi-|  that 


inf  P„(CS|S2)  -  inf  P„(Xk.,*,.kIP0  - 


=  P(Vi  >  -  d). 


We  need  the  constant  d  such  that  P{M^  >  -  d)  >  P*  holds  for  all  i, 

1  £  i  <_  k.  By  Theorem  3.1  we  have  d  =  dj]j|.  It  also  follows  that 
if  and  Sg  are  the  selected  subsets  associated  with  selection  proce¬ 
dures  «i  and  fig,  respectively,  then  S1  c  Sg.  Thus  61  is  better  than  6^,. 


Definition  3.4.  The  procedure  fi^  is  defined  as  follows:  Let  Xj-max(X^,. 
Step  1.  Select  7^.,  i  1  and  stop,  if 

*'-’*•*'  f 

otherwise  reject  tt^  and  go  to  Step  2. 

Select  tt^  ,  i  >_  2  and  stop,  if 


Step  2. 
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otherwise  reject  n0  and  go  to  Step  3. 


Step  k-1.  Select  ^ .  i  1  k  -  1  and  stop,  if 

Xk.1  1  >'0  -  ^ 

otherwise  reject  and  go  to  Step  k 
Step  k.  Select  ^  and  stop,  if 

*k  ~  V°  dk  ^  ’ 

otherwise  reject 

Hered.'s  are  the  smallest  values  such  that  63  satisfies  the  P*-condition. 

Evaluation  of  d.  's 

For  any  i ,  1  £  i  <  k, 

inf  P  (CS|6_)  = 
u€Oi  £  * 

k-i+1  - 
-  Pu*(  0  {X 
y  j=i  J 

=  p(zk-i+l  -  ‘  dk-i+l* 

whenever  y*  =  (uq>-°°,. . . ,-®,wq,. . .  ,uq)  €  . 

i 

Since  is  N(0 , 1 ) ,  it  implies  dk_i+1  =  d  for  all  i,  and 
d  «  $_1(P*). 

Hence,  we  have  the  following  theorem: 


k-i+1  .. 

inf  P  (  U  (X.  >  pn  -  d.  -— }) 
i!-  j=l  J  0  J  v^T 


-  i 


-rl  .  \ 
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Theorem  3.5.  Selection  procedure  satisfies  the  P*-condition  with 
di  =  d,  i  =  which  do  not  depend  on  i.  Hence  the  procedure  is 

not  changed  if  the  statistics  are  replaced  by  ,  the  sample  mean  of 
population  tt .  for  i  =  l,...,k. 

The  following  procedure  6^  was  given  by  Gupta  and  Sobel  (1958), 
without  assuming  any  ordering  prior: 

Definition  3.5.  The  selection  procedure  6^  is  defined  as  follows: 

o  •  Select  it.  if  and  only  if  X.  >  -  d  —  i  = 

4  1  1  ~  0  /n 

where  d  is  the  smallest  constant  such  that  <5^  satisfies  the  P*-condition. 

It  was  shown  that  the  value  d  is  determined  by  the  equation 

1  1 

*(-  d)  =  1  -  P*^  i.e.  d  =  $"'(P*k). 

(?) 

3.5.  Some  Proposed  Selection  Procedures  Si  i  =  1,  2,  3,  4 
When  y q  is  Unknown 

p 

Case  II.  Uq  unknown,  connion  a  known,  common  sample  size  n. 


Definition  3.6.  We  define  a  selection  procedure  6^'  by  replacing 
the  inequalities 


Xi:k  -  y0  "  di:k 


(Di  i  =  )  k 


in  procedure  6-j  (Definition  3.1)  with 

X.  .k  >_  XQ  -  d|?k  •— ,  i  =  1 , . . .  ,k,  respectively. 


Here  XQ  =  J  XQl./n,  d|fj',  i  =  l,...,k  are  the  smallest  constants  such 
that  the  selection  procedure  6^  satisfies  the  P*-condition. 

Similar  to  the  Case  I,  we  have  the  following  theorem: 
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Theorem  3.6.  For  any  1 ,  1  <  1  <  k,  d^]+1.k  Is  determined  by  the 
equation 

/  P(V1  1  t  -  d^+i:k)d,(t)  -  P*.  (3.15) 

«oo 

It  is  easy  to  see  that  d^|+1,k  =  dj?j  and  it  Is  Increasing  In  i.  The 

1  (2) 
lowing  theorem  gives  us  an  identical  form  of  the  selection  procedure  '. 

Theorem  3,7,  The  selection  procedure  is  not  changed  if 

A 

the  statistics  Xi-k,  i  =  l,...,k,  are  replaced  by  X1>k,  i  *  l,...,k, 
respectively. 

Proof.  The  proof  is  straightforward  and  hence  it  is  omitted. 

The  values  dj?j,  i  =  l,...,k  are  tabulated  in  Table  II  for 
k  =  1  (1)  6,  8,  10,  -  and  P*  =  .99,  .975,  .95,  .925,  .90. 

Similar  to  the  Case  I,  we  propose  a  selection  procedure  as 
follows: 

Definition  3.7.  We  define  a  selection  procedure  6 ^  by 

<$2^:  Select  ttj  If  and  only  if  X^.k  >_  Xq  -  d  —  i  =  l,...,k 

/n 

where  d  is  the  smallest  value  such  that  6^  satisfies  the  P*-condition. 
Then,  similar  to  procedure  we  have  d  =  d^. 

Next,  we  define  a  selection  orocedure  which  is  similar  to  <$3. 
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Definition  3.8.  The  selection  procedure  <5^  is  defined  by  replacing 

>.  -  d.  ~  In  $3  (Definition  3.4)  by  >  Xq  -  d|  —  ,  1  =  1 ,.. .  >k 

</n  /n 

(2) 

where  d|,...,d^  are  the  smallest  values  such  that  «3  satisfies  the 
P*-condition. 

Similar  to  Theorem  3.5  we  have: 

Theorem  3.8.  The  selection  procedure  6^  satisfies  the  P*-condition 
with  di  =  d,  i  =  l,...,k  where  d  is  determined  by  the  equation 

rr> 

/  *(d-t)d*(t)  =  P*.  (3.16) 

—  00 

And  6^  is  not  changed  if  the  statistics  is  replaced  by  , 
the  sample  mean  of  population  n..  for  i  s  l,...,k. 

(2) 

The  following  selection  procedure  6^  '  was  proposed  by  Gupta  and 
Sobel  (1958): 

Definition  3.9.  The  selection  procedure  6^  is  defined  by 

Select  r,  1*  and  only  if  >  Xq  -  d  — -  i  =  1,...,k 

^7 

where  d  is  determined  by  the  following  equation. 

/  n  0(u  f^r  +  d)]$(u)du  s  P*. 

—  i«l  /  "0 


(3.17) 
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For  the  special  case  *  n  (1  *  0,  k) 

/  * k(t+d) 4»(t)dt  *  P*. 


(3.18) 


Under  the  normal  distribution  N(0,1),  the  tables  of  d-values  sat¬ 
isfying  the  Equation  (3.18)  for  several  values  of  P*  are  given  In 
Bechhofer  (1954)  for  k  =  1  (1)  10  and  in  Gupta  (1956)  for  k  *  1  (1)  50. 

3.6.  Some  Proposed  Selection  Procedures  «j3\  1  *  1,  2,  3,  4  for  the  Normal 

2 

Means  Problem  When  Common  Variance  a  Is  Unknown 

2 

Case  III.  Uq  known,  common  variance  o  unknown,  n^  *  n  >  1. 


Definition  3.10.  We  define  the  selection  procedure  fij3'  by  replacing 
the  inequalities 

Xt:ki',0  ■  di!k^  1  *  1 . k 

in  procedure  6^  (Definition  3.1)  by 
A  (3)  S 

Xi:k  -  w0  '  di  :k  T-  *  =  respectively, 

where  d^'s  are  the  smallest  values  such  that  satisfies  the 
P*-cond1tion  ;  S  denotes  the  pooled  estimator  of  o  based  on 


v  *  k(n-l),  that  is 


k  n 


s  -  l  I  (Xn  -  XJVv. 
1=1  j=l  tJ  1 


(3.19) 
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Note  that  -jjj-  has  the  chi-square  distribution  xv  with  v  degrees  of 

0 

freedom.  The  following  theorem  then  follows: 


Theorem  3.  9. 
is 

or 


(3) 

The  equation  which  determines  the  constant  dk_.j+^.k 


P<V1  i  -  "k-Hl  :k  f>  ■  P* 


(3.20  ) 


l0  P(Vi  dk-V*l :k  ■  p* 


(  3.21  ) 


where  qv(y)  is  the  density  of  . 


We  can  rewrite  Formula  (3.21)  as 

L  p<*,  i  - /!><(«  ■ p* 


or 


-  -1 
o  1 


//(*,  i  -  :k  /?)  <“  * p*- 


(3.22) 


‘2  T" 
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Remark 


3. 2.  The  values  of  dj^]+1.k,  1  ■  l,...,k  depend  on  \>  ■  k(n-l);  also 


d(3)  y 

°k-1+l:k  r  al:T 


By  using  Rablnowltz  and  Weiss  table  (1959)  (w1thNa24  and  n  of  their  table 

k-1+1  :k*  1‘1 . k* 


equal  to  0)  we  have  evaluated  and  tabulated  the  values  of  d^ 


In  Table  III,  for  k  -  2  (1)  6,  P*  «  .99,  .975,  .95,  .925,  .90, 
with  common  sample  size  n  =  3,  5,  9,  and  21. 

For  k  >_  6  and  n  >  21  ,  l.e.  v  >  120  we  can  reasonably  well  approximate 


,(3) 

3k-1+l  :k 


by  di; 


(1) 


1:1 


Definition  3.  11.  We  define  the  selection  procedure  <5^  by 

6^:  Select  if  and  only  If  ,k  >_  yQ  -  d^  —  1  =  1 ,. . .  ,k 

Jr\ 

where  S  is  defined  as  In  procedure  and  Is  the  smallest 
constant  such  that  satisfies  the  P*-cond1t1on. 

As  before,  it  can  be  shown  that  d^  *  d|^. 

(3) 

Remark  3.3.  In  Case  III  the  selection  procedure  6]  will  not  be 

changed  if  we  replace  the  isotonic  statistics  X^.k  by  X^,  respectively. 

(31 

But  this  is  not  necessarily  true  for  selection  procedure  . 


Definition  3,12.  The  selection  procedure  6^  Is  defined  to  have  the 
same  form  as  procedure  6^  except  that  the  inequality  defined  In  the 

(2) 

ith  step  of  procedure  6^ '  is  replaced  by 

X.j  >  Uq  -  d  —  for  1  =  l,...,k. 


The  proof  of  the  following  theorem  uses  the  same  arguments 
as  that  In  Case  I,  hence  it  is  omitted. 


Theorem  3.10.  The  equation  which  determines  the  constant  d  of 
selection  procedure  6^  is 

/  $(yd)qv(y)dy  *  P*.  (3.23) 

Gupta  and  Sobel  (1958)  gave  a  selection  procedure  6^  '  in  this 
case.  It  is  as  follows: 

Select  tt,.  if  and  only  if  >_  Uq  -  D  — —  i  =  l,...,k 

and  the  equation  which  determines  D  Is 

^  *k(yO)qv(y)dy  =  p*.  (3.24) 

where  v  *  l  (n,-l). 
i*l  1 

3.7.  Some  Proposed  Selection  Procedures  i  3  1,  2,  3,  4  for  the  Normal 

2 

Means  Problem  When  Both  Control  and  Common  Variance  a  are  Unknown 

2 

Case  IV.  pq  unknown,  common  variance  a  unknown  and  common  sample 
size  n. 

(3) 

Here  we  replace  jjq  in  each  selection  procedure  by  XQ,  1  <  j  <  4, 
and  get  procedures  1  <  j  <  4,  respectively.  The  constants  ;|<> 

i  =  l,...,k,  of  procedure  are  determined  by 


24 


y 


f  f  P(*.  i  a  -  4-i+l  .|(  /|)d*(u)dx^(t)  =  P*.  (  3.  25) 

0  -00 

The  constant  d  of  procedure 


The  constants  d  of  procedures  6^  and  6^  are  determined  by 

f  /V(u^d)d,(u)dxf(t)  *  P*  (3.26) 

0  -»  v 

with  r  =  1  and  k,  respectively,  and  their  values  for  selected  values 
of  P*,  k  and  v  are  given  in  Gupta  and  Sobel  (1957)  and  Dunnett  (1955). 

3.8.  Properties  of  the  Selection  Procedures 
Under  simple  ordering  prior,  it  is  natural  to  require  that  an  ideal 
selection  procedure  is.  isotonic  as  defined  below: 

Definition  3.13.  A  selection  procedure  6  is  isotonic  ii  it 
selects  tt^  with  parameter  ,  and  if  <  pj,  then  it  also  selects  ifj. 
Procedure  6  is  weak  isotonic  or  monotone  if 

F(tt.  is  selectedjs)  £  P(ir.  is  selected|s)  whenever  y.  <  y.. 

*  J  '  J 

It  is  easy  to  see  that  any  isotonic  selection  procedure 
is  weak  isotonic,  but  the  converse  is  not  true. 

Now,  let  6  P  ^  s  5.,  i  =  1,  2,  3,  4. 

Theorem  3. 11  .  The  selection  procedures  <sj^,  and  6^  are 
isotonic  and  procedure  ^  is.  monotone,  for  i  =  1,  2,  3,  4. 


j.  j  .  ri  it  — * 


m 

ZL 


'''  ''"'^pppppi^p 
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Proof.  The  proof  follows  immediately  from  the  definitions  of  the  pro¬ 
cedures. 

Given  observations  1  =  x  =  (Xq,...,x^)  where  x..  is  the  sample  mean 
of  population  ,  i  =  l,...,k,  and  xQ  -  pQ  if  pQ  is  known,  otherwise 
Xg  is  the  sample  mean  of  population  ttq.  Let 

4^(x,  $)  =  PCtt^  included  in  the  selected  subset|X  =  <$) 

for  i  =  1 , . . . ,k. 

Definition  3.14.  A  selection  procedure  S  is  called  translation- 

k+1 

invariant  if  for  any  x  t  F  ,  c£R 
^(Xg  +  c,  x1  +  c,...,xk  +  c;  6)  =  ^(xg, - x.  ;  s),  i  =  1 ,. . .  ,k. 

Theorem  3.12.  The  selection  procedures  6^,  6^  and  <s^ 

are  translation-invariant  for  i  =  1,  2,  3,  4. 

Proof.  Proof  is  straightforward  and  hence  omitted. 

Expected  Number  (Size)of  Bad  Populations  in  the  Selected  Subset 

Suppose  the  control  pg  is  known  and  we  have  common  sample  size  n 

2 

and  common  known  variance  a  ;  without  loss  of  generality,  we  assume 
that  pg  =  0  and  o/vW  =  1.  Let  E( S 1 ] 6 )  denote  the  expected  number  of 
bad  populations  in  the  selected  subset  in  using  the  selection  procedure 
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On  the  other  hand,  for  procedure  <$£ 

sup  E{S'1«2)  =  I  p(  S  {Z,.i  1-  d!!jb)  .  (3.28) 

p€Pk_j  ■  r=l  1=1  , 

From  (3.28)  we  see  that  the  supremum  for  1S  increasing  in  j  and 
is  greater  than  or  equal  to  the  supremum  for  61  given  in  (3.27),  since 


d(l)  =  d(l)  <  dO) 

°t:k  °l:k-t+l  -  °l:k' 


Therefore,  we  have  the  following  theorem  (see  also  the  remark  just  before 
Def,  3.4). 

Theorem  3.13.  For  any  i »  0  £  i  £  k 

sup  E(S '  | )  >  sup  •KS*  16,), 

iLeni  6  p€ni 


sup  E( S  * | )  =  sup  E(S' j  6p ) . 


Theorem  3.14.  In  Section  3.1,  Case  I,  for  any  j,  0  <_  j  <_  k 

sup  E(S'|s,)  -  j  -  q(l-qJ)/P*  (3.29) 

Ji&k-j 


where  q  =  1  -  P*. 


"T . „ 


. . :  ^  ■  ”•  -  v"?\\  '^T^:r»fr:Ti  v.  A>\tr  '*■*■  "vv n  JilJIfR'UJl 


Proof. 


sup  E(S' |d,) 
— ^k-j 


=  sup  J  P  (select  ». |j,) 


=  sup  z  P  (  ma*  X  >_  -  d) 
jj€o  i=1  U-  l<r<1  r 

=  SO-  n  F(-d))  " 

i=l  r-1  N  ■ 

=  j  -  j  q1 
1=1 


=  j  -  q(l-qj)/P* 


where  q  =  (1-P*). 


Theorem  3.15.  sup  E(S'!<$,)  is  increasing  in  j,  hence 


p«nk-j 


sup  E(S'|«,)  -  sup  E(S' |«,)  =  k  -  ql .-Uk!/P*.  (3.30) 


^k-j 


Proof.  Since 


(j+1) 


(j-  {  q1)  -  l-qj+1  >  0. 

i=l 


In  Case  I  , 


Gupta  (1965)  showed  that 


sup  E(S ' |fi4) 


1 

kP*k. 


(3.31) 


1 


tUhV'^lh  I'iiukkr  .lil.fj 


1— S  r  C Jit 


%  p'lluea  I 


— ”7T 


it.  BAtOi  it  uki  jglMi&s 
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Let  us  define  the  event  A1  =  (Z1;k  >  -  1  *  1 _ _ ,k;  then 

we  have 

J  j 

Lemma  3.2.  P(  u  A,  n  A.,,)  >  P(  u  AJ  P*  for  all  j,  1  <  j  <  k-1, 
1=1  1  J  1  1=1  1  -  _ 


Proof:  P(  u  A,  flAul) 
i=l  1  J  1 

=  P(1Vi1:ji-d‘!k,n  Vl> 

>  P(  ut  A, )P(AJ+1 ) 

-  P(  U  A.)  P*. 

1=1  1 


The  above  inequality  is  a  result  of  the  fact 

A  . 

Ai  c  {^i:j  -  "  d1:k}  for  a11  1  "  j  =  1 . k-1. 

Theorem  3.16.  For  all  k  2,  sup  E( S ’  | <5 )  <  sup  E(S ' 1 6->) . 

f‘0  fi0 

Proof:  To  prove  the  theorem  it  is  sufficient  to  show  that  for  all 

j  i 

given  k  >  2,  P(^j  A^)  <  1  -  (1-P*)J  for  all  j  and  strictly  inequality 
holds  for  some  j,  1  <_  j  ^  k. 

It  holds  for  j  =  1,  since  P(A,)  =  P*.  Suppose  P(  u  A. )  <  1-(1-P*)^ 

1  1-1  '  ~ 
is  true  for  some  j,  1  <  j  <  k-1,  then 
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j+1  J  j 

P(  U  Ai )  =  P(  y  A.)  +  P*  -  P(  U  A.  n  A,  ,) 
i«l  1  1=1  1  1=1  1  J+l 

j  j 

<  P(  u  A.)  +  P*  -  P(  U  A. )P* 

1=1  1  1=1  1  | 

<  P*  +  (1-P*)(l-(1-P*)j)  j 

i 

=  1  -  0-p*)j+1.  j 

j 

Hence  by  induction  principle,  the  proof  is  finished. 

i 

This  theorem  tells  us  that  procedure  6^  is  better  than  63  in  j 

the  sense  that  in  Uq  it  tends  to  select  smaller  number  of  bad  populations, 
however,  procedure  5-j  is  not  uniformly  better  than  6j.  In  some  cases 
(see  Section  5),  is  slightly  better  than 

When  the  ordering  prior  among  the  unknown  parameters  is  unknown, 
we  can  use  the  selection  procedure  of  Gupta  and  Sobel  (1958)  or  use  j 

I 

the  ordering  of  the  sample  means  as  the  ordering  of  unknown  parameters  5 

] 

and  apply  the  selection  procedure  which  is  originally  used  under  ordering  prior.  ' 

In  the  normal  case  with  the  latter  approach,  the  substitution  implies  that  the 

J 

isotonic  regression  of  the  sample  means  turns  to  the  usual  ordered  j 

sample  means,  and  that  the  selection  procedures  6^ i  =  1 ,  2,  3,  4,  j 

are  of  the  same  type  as  '  (i  =  1,2, 3, 4),  respectively,  and  the  selection  j 

procedures  <5^,  j  =  1,3,  i  =  1,2, 3, 4  are  of  the  same  form  as  65^,  i=l  ,2,3,4,  j 
respectively,  which  are  equivalent  to  the  procedures  proposed  by  Naik  j 

(1975)  and  Brostrom  (1977),  independently  (see  also  Holm  (1979)).  1 

A-  Selection  Rules  for  the  Location  Parameter  Under  Partial 
Ordering  Prior  Assumption 

Assume  that  we  have  only  a  partial  ordering  prior  of  k  unknown 
location  parameters,  that  is  the  parameter  space  j 


L 

ft'a  {o|e€F  and  there  Is  a  partial  order  relation  "<"  among  e^'s) 
Our  approach  Is  to  partition  the  set  (e1 . . ,ek}  Into  several  sub- 

SL 


sets,  say  B 


O’* 


,Bf,  so  that  B 


n  Bj  * 


0,  If  1  t  j, 


U  B. 
j=l  j 


and  for  each  Bj  {j  =  1, there  is  a  simple  order  on  It  and  there 
is  no  order  relation  among  the  elements  of  subset  BQ. 

Let  *  |B.j  | ,  the  number  of  elements  contained  in  B. ,  1  =  0, 
so  we  have 


i 

E  b  =  k. 

1*0  1 

If  we  denote  the  new  Induced  partial  order  by  then  we  have 

a  parameter  space  ft"  =>  ft'.  We  use  an  example  to  Illustrate  how  to  find 
an  Induced  partial  order. 


Example.  Suppose  k  =  8,  and  we  have  a  partial  ordering  prior  <  e^, 

91  -  °8’  91  -  e2  -  93  -  94’  anc*  e2  -  96  -  e7'  We  use  a  "tree"  to 
represent  this  partial  ordering  as  in  Figure  1. 


Figure  1.  Original  partial  ordering 
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Then  we  have  an  Induced  partial  ordering  <  8?  <  e3  <  04,  eg  <  07 
as  In  Figure  2. 


e4r 


03  * 


V 


Figure  2.  Induced  partial  ordering. 

And 

B0  =  {05’  08} 

Bi  =  ©2»  e4^ 

B g  “  1 6g  *  ^7 1 • 

It  is  clear  that  the  induced  partial  order  is  not  unique,  for 
example,  we  can  partition  {o^,...,0gi  into  three  other  subsets 
Bq,  Bj ,  B£  where 

B0  =  {e5*  e8} 

Bj  =  {0-j ,  ©2 ,  »  67) 

B2  =  {03’  V' 

For  the  location  parameter  case,  a  selection  procedure  sp  can  be 
defined  as  follows: 

Definition  4,i.  We  define  a  selection  procedure  6P  as  follows: 

Suppose  Bq,...,B^  are  induced  subsets  and  that  for  each  subset 

B.,  j  =  l,...,t  there  is  a  simple  order  on  it.  We  choose  a  proper 
J 


where  p.g  is  the  parameter  space  associated  with  the  subset  B^. 

5.  Comparisons  of  the  Performance  of  Basic  Rules  for  the  Normal  Means  Problem 

In  this  section  we  describe  results  of  a  Monte  Carlo  study  to  compare 

the  performance  of  selection  procedures  5j,  62,  6y  and  6^..  Suppose  we  j 

2  j 

have  k  independent  populations,  each  population  with  distribution  N(uit  o  ),  j 

2 

with  common  known  variance  u  and  coimion  sample  size  n.  Assume  that 

the  mean  uq  of  the  control  is  known;  without  loss  of  generality  we  * 

assume  that  Uq  =  0  and  o/v'n  *  1.  • 
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In  the  simulation  study,  we  used  Rubin  and  Hinkle's  RVP-Random  Variable 
Package,  Purdue  University  Computing  Center,  to  generate  random 
numbers.  For  each  k,  we  generated  one  random  number  (variable)  for 
each  population,  then  applied  each  selection  procedure  separately 
and  repeated  it  ten  thousand  times;  we  used  the  relative  frequencies 
as  an  approximation  of  the  exact  values  of  the  associated  performance 
characteristics  for  each  procedure.  In  Table  IV  we  use  the  following 
notations: 

u  =  (u^,...,uk),  ik  is  the  parameter  of  population  n.. 

PS  =  P(CS) 

PI  *  P(correctly  rejecting  ail  bad  populations) 

PC  =  P(correct  classification  of  all  popuiuL.on) 

where  the  correct  classification  means  that  we  select  all  good 

populations  and  reject  all  bad  populations. 

El  =  Expected  number  (size)  of  bad  populations  contained  in 
the  selected  subset. 


Ed 


l  <i 


Vy0 


P  (ir.j  is  selected) 


ES  =  Expected  size  of  the  selected  subset. 


Table  IV. 1  consists  of  four  parts,  namely,  the  four  values  of 
k  =  2, 3,4,5,  for  each  value  of  k  we  assume  that  we  have  two  bad  populations. 
In  this  case  based  on  the  performance  characteristics  PI,  PC,  El  or  EJ,  we 
found  the  performance  ordering  as  follows: 

^i  >■  ^2  >  ^3  v  ^4* 

where  <$i  >  6^  means  that  <$.j  is  better  than 
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In  Table  IV. 2  we  assume  that  we  have  three  bad  populations  for  k  =  3, 
and  that  both  populations  are  bad  for  k  *  2,  this  table  indicates  the  same 
trend  as  Table  IV. 1,  i.e.  5-|  >.  •>  v  6^.  If  k  is  increased  by  adding 
strictly  good  (parameter  strictly  larger  than  control)  populations,  then 

ft  A 

A  * 

E I { 6 ,  1  =  1,2  does  not  increase.  This  Is  because  X^.k  >  X^,k+^  a.s. 

1  <  1  <  k. 

In  Table  IV. 3  we  assume  that  for  each  k,  k  =  2, 3, 4, 5  that  every  population 
is  bad.  Based  on  the  quantities  PI,  PC,  El  and  EJ,  we  find  that  the  perfor¬ 
mance  is  as  follows: 

4  -j  v  62s-  63  >  4^. 

This  is  the  same  result  as  before. 

Table  IV. 4  has  the  same  structure  as  before,  but  for  each  value  of  k, 
k  =  2, 3, 4, 5,  we  assume  that  the  first  population  is  the  one  and  only  one 
bad  population  with  parameter  -1  which  is  less  than  the  control  =  0. 

A  glance  at  the  table  indicates  that  the  performance,  based  on  the  charac¬ 
teristics  PI,  PC,  El  and  ES,  can  roughly  be  ordered  as  follows: 

63  ■>  (Sg*  ]  v  4^. 

i.e.  procedure  63  is  the  best  and  is  slightly  better  than  6^  ancl  6i>  52 
and  6-j  are  very  close  and  both  are  better  than  6^.  As  the  number  of 
populations  k  increases  from  two  to  five  and  the  three  additional  populations 
are  good  populations  with  parameter  1,  2,  and  3,  respectively,  we  find 
that  El (5i ,  k  «  5)  -  E I ( 6 i ,  k  =  2),  i  =  1,2, 3, 4,  is  0.0124,  0.0124,  0.0031, 
0.121,  respectively.  This  means  that  when  k  increases  and  the  additional 
populations  are  good,  then  procedure  6^  is  the  most  sensitive  procedure  with 
k  and  thus  not  good  in  terms  of  El.  6^  seems  to  perform  better  in  terms  of 


El  while  ^  and  62  are  about  the  same. 
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In  Table  IV. 5  we  assume  that  the  ordering  prior  of  unknown  parameter 
Is  Incorrect;  i.e.  the  true  configuration  (-2,  -1,0, 1,2)  is  replaced  by 
(-1 ,-2,1 ,0,2) .  The  simulation  results  indicate  that,  based  on  PI,  PC, 

El  and  EJ  we  have  performance  6-j  >  <$2  >  [63 . <S^] .  Thus  here  again  61  is 
the  best.  If  we  compare  Table  IV. 5  with  Table  IV. 1,  we  see  that  6^  does 
not  change  (the  small  differences  are  because  of  random  fluctuations), 
EI(63)  and  EJ{<53)  increase  quite  appreciably. 

From  these  five  tables,  it  appears  that,  in  general,  the  overall 
performance  of  these  procedures  is  <5^  6^  >  63  >  <5^,  if  the  ordering 
prior  is  correct.  If  there  is  no  information  regarding  the  prior  ordering, 
then  or  <Sg  seem  to  be  an  appropriate  procedure  to  use, 
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TABLE  I 

Table  of  d^}},.  values  (satisfying  (3.5)  and  (3.14))  necessary  to 
carry  out  the  procedure  for  the  normal  means  problem  under  the  simple 
ordering  prior. 


\d(l> 

N\l:k 

k 

* 

P 

.99 

.975 

.95 

.925 

.90 

1 

2.3264 

1.9600 

1.6449 

1.4395 

1.2816 

2 

2.3337 

1.9775 

1.6780 

1.4872 

1.3430 

3 

2.3339 

1.9787 

1.6817 

1.4942 

1.3538 

4 

2.3339 

1.9787 

1.6823 

1.4956 

1.3563 

5 

2.3339 

1.9787 

1.6824 

1.4960 

1.3571 

6 

2.3339 

1.9787 

1.6824 

1.4960 

1.3o73 

00 

2.3340 

1.9787 

1.6824 

_ 

1.4960 

1.3574 

TABLE  II 


(2) 

Table  of  values  (satisfying  (3.15))  necessary  to  carry  out  the 

(2) 

procedure  6v-j'  for  the  normal  means  problem  under  simple  ordering  prior. 


\  d(2> 
n|i  :k 

k 

★ 

P 

.99 

...  . 

.975 

.95 

.925 

.90 

1 

3.2886 

2.7711 

2.3258 

2.0355 

1.8122 

2 

3.3449 

2.8494 

2.4267 

2.1530 

1.9434 

3 

3.3605 

2.8730 

2.4589 

2.1917 

1.9874 

4 

3.3673 

2.8840 

2.4723 

2.2105 

2.0091 

5 

3.3711 

2.8901 

2.4832 

2.2215 

2.0219 

6 

3.3734 

2.8941 

2.4890 

2.2286 

2.0303 

8 

3.3761 

2.8988 

2.4960 

2.2375 

2.0406 

10 

3.3776 

2.9014 

2.5000 

2.2426 

2.0440 

3.3787 

2.9032 

2.5021 

2.2448 

2.0487 
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TABLE  IV. 2 


|  Simulation  results  for  the  comparative  performance  of  various  selection 
|  ■  procedures  for  the  normal  means  problem  (notation  explained  in  Section 
I  5)  under  simple  ordering  prior. 
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TABLE  IV. 3 


Simulation  results  for  the  comparative  performance  of  various  selection 
procedures  for  the  normal  means  problem  (notation  explained  in  Section 
5)  under  simple  ordering  prior. 
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TABLE  IV. 4 


Simulation  results  for  the  comparative  performance  of  various  selection 
procedures  for  the  normal  means  problem  (notation  explained  In  Section 
5)  under  simple  ordering  prior. 
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Simulation  results  for  the  comparative  performance  of  various  selection 
procedures  for  the  normal  means  problem  (notation  explained  in  Section 
5)  under  simple  ordering  prior. 
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